Particle number renormalization in almost half filled Mott Hubbard superconductors 
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The effects of the Gutzwiller projection on a BCS wave function with varying particle number are 
considered. We show that a fugacity factor has to be introduced in these wave functions when they 
are Gutzwiller projected, and derive an expression for this factor within the Gutzwiller approxima- 
tion. We examine the effects of the projection operator on BCS wave functions by calculating the 
average number of particles before and after projection. We also calculate particle number fluctua- 
tions in a projected BCS state. Finally, we point out the differences between projecting BCS wave 
functions in the canonical and grand canonical schemes, and discuss the relevance of our results for 
variational Monte Carlo studies. 
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PACS numbers: 74.20.-z, 74.20.Mn 



I. INTRODUCTION 



Recently, Anderson has underscored the importance of 
fugacity in wave functions that do not conserve particle 
number 1 . Following an earlier paper by Laughlin-, An- 
derson argued that a fugacity factor should be included 
in variational wave functions of the form, 



P I* 



BCS/ 



P 



n 



Uk 



Wfe4t ct 



(1) 



Here, P — Yii(^~ n it n il) is a (Gutzwiller) projection op- 
erator which excludes double occupancies at sites i (Ref. 
3), and |\&bcs), a BCS wave function. Projected wave 
functions of this form were originally proposed to de- 
scribe the phase diagram of doped Mott Hubbard insula- 
tors such as the high temperature superconductors^*^. 
Detailed variational Monte Carlo (VMC) studies have 
been carried out recently using projected d-wave BCS 
states as variational wave functions for the two di- 
mensional Hubbard model, after a suitable canonical 
transformation 



Despite their simple form, projected wave functions 
exhibit nontrivial properties because the projection op- 
erator acts on a quantum many body state. The ac- 
tion of the projection operator (in reducing the allowed 
states in the Hilbert space) concomitant with the corre- 
lations of the quantum state being projected, leads to a 
variety of physical phenomena^. Other interesting ef- 
fects include non-trivial matrix-element renormalization 
near half- filling 11 and the occurrence of superconductiv- 
ity near an (antiferromagnetic) Mott insulator—. 

Approximate analytical calculations with wave func- 
tions such as Eq. Q can be done using a renormalization 
scheme based on the Gutzwiller approximation. Within 
this approximation, the effects of projection on the state 
|\&o) are approximated by a classical statistical weight 
factor multiplying the quantum resuMpi&i^. Thus, for ex- 



ample, 
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where O is any operator, and g, the so called Gutzwiller 
factor. For example, the Gutzwiller approximation for 
the kinetic energy operator c\cj + c|cj and the superex- 

change interaction between sites i and j, Si- Sj yields the 
Gutzwiller factors, 
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where n is the density of electrons. In deriving these 
renormalization factors, one considers the number of 
terms that contribute to (*|6|*) and to (* |O|* ) re- 
spectively. The ratio of these two contributions is the 
renormalization factor. 

The renormalization factors are functions of the local 
charge density. This is a well defined quantity when one 
considers for example, a projected Fermi liquid state, 
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But suppose instead, we consider wave functions such 
as the BCS state in Eq. (JJ, where the particle number 
fluctuates. In this case, it is not clear what the local 
charge density in Eq. @ should be. It may be argued 
that the correct n is set by the average particle number 
N. But then, projecting a BCS state changes the average 
particle number; i.e., the average number of electrons in 
I^bcs) does not equal that in P^bcs)- Clearly, we need 
a scheme to keep track of this effect. 

Note that this problem can be avoided completely, as 
is done in most variational Monte Carlo (VMC) stud- 
ies. Here, the particle number is fixed (one works in a 
canonical ensemble), and Eq. replaced by, 
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The operator P/v fixes the particle number, and the issue 
of projection changing the mean particle number does not 
arise^. However, there are also other VMC studies with 
wave functions that do not have fixed particle number—. 
Moreover, we are often interested in carrying out analyt- 
ical approximations in the spirit of Eq. J5J. Since such 
manipulations are easier done with BCS wave functions 
(where the particle number is not fixed), it is desirable 
to understand the effects of the projection operator on 
this class of wave functions. In this paper, we present 
analytical and numerical considerations of this problem. 
In doing so, we clarify the notion of fugacity introduced 
by Andersoni. We also discuss the relevance of this ap- 
proach for the Gutzwiller approximation in the grand 
canonical scheme and the corresponding VMC studies^. 



II. THE FUGACITY FACTOR 

Consider the projected BCS wave function, Eq. QJ. 
It is clear that the projection operator P changes the 
average number, viz., 
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The effect of the projection operator can be seen most 
clearly by examining the particle number distribution in 
the unprojected and projected Hilbert spaces. Towards 
this end, let us write the average numbers, N( QS> (N) in 
the unprojected (projected) Hilbert space 
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are the particle number distributions in the unprojected 
and projected BCS wave functions respectively; Pn is an 
operator which projects onto terms with particle number 
N. The particle number distributions before and after 
projection may be related by 
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Eq. (|SJ) constitutes the Gutzwiller approximation for the 
projection operator Pn with the corresponding renormal- 
ization factor, gN] C is an irrelevant constant (the ratio 
of the normalization of the unprojected and projected 
wave functions), which does not depend on N. Follow- 
ing Gutzwiller, we estimate qn by combinatorial means, 
as being equal to the ratio of the relative sizes of the 
projected and unprojected Hilbert spaces. Then, 
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where L is the number of lattice sites and iVf (iVj) is 
the number of up (down)-spins. Since in a BCS wave 



function, TV 



N 



N/2, N being the total number of 



particles, the expression for gN can be simplified to 
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Hence, if we were to impose the condition that the av- 
erage particle number before and after projection be iden- 
tical, a factor g^ 1 has to be included in Eq. Q. Then, 
from Eq. Q and Eq. (JSJ, we obtain the particle number 
after projection N ncw , 
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which is the desired result. 

Now, let us show how this procedure can be imple- 
mented for the wave function I^bcs)- Since the BCS 
wave function is a linear superposition of states with par- 
ticle number . . . , N — 2, N, N + 2, . . ., we consider the ef- 
fect of projection on two states whose particle numbers 
differ by 2. Then, the ratio, 



9N+2 
9N 



L-N 
L-N/2 



(12) 



in the thermodynamic limit. Eq. 1121 shows that the pro- 
jection operator acts unequally on the N and 7V+2 parti- 
cle states; the renormalization of the weight of the N + 2 
particle states gN+2, is f 2 times the weight of the N par- 
ticle states, gN- This effect can be rectified as in Eq. <|11|) 
by multiplying every Cooper pair c^clj.| by a factor j 
in the BCS wave function. It produces the desired result, 
viz., the projected and unprojected BCS wave functions 
have the same average particle number. 

Alternatively (following Anderson), we can multiply 
every empty state by the factor / and write, 
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Then again by construction, the fugacity factor / in 
Eq. (|13|) ensures that the projected wave function 
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P l^jfcg) and the unprojected wave function |\I/bcs) have where, 
the same particle number. The denominator in Eq. I|13|) 
is the new normalization factor. 

The following points are in order: (a) the fugacity 
factor / in Eq. (|12[1 depends on the variable particle 
number N . However, since the particle number of the 
BCS wave function is sharply peaked within the range, 
N(o) _ ViVW and N<® + VW), we will assume that the 
fugacity factor / = /(A^ )) in the thermodynamic limit; 



g;, where g t is 



(b) in this limit, Eq. (jT2|) reduces to f 2 
the Gutzwiller factor defined in Eq. ©. Then, Eq. fT3|) 
reduces to 
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which is the wave function proposed by Andersoni; (c) 
the fugacity factor ensures that projection affects the N 
and N + 2 particle states of the BCS wave function in the 
same way. In principle, such a factor could depend on the 
k- value, but in this paper, we will treat it as a mere com- 
binatorial device; (d) the combinatorial argument fails at 
half filling when L = N<-°\ 



III. PARTICLE NUMBER RENORMALIZATION 
IN PROJECTED BCS WAVE FUNCTIONS 



In the previous section, we showed that the inclusion 
of the fugacity factor is necessary for the average parti- 
cle number in a BCS wave function to be unchanged by 
projection. Alternatively, one might ask what is the ef- 
fect of the projection operator on a BCS wave function; 
viz., if projection changes the mean particle number of 
a BCS state, how are the particle numbers before and 
after projection related? In this section, we will use the 
fugacity factor to answer this question. In particular, 
we will show how particle density after projection can 
be determined self consistently by including the fugacity 
factor in the projected BCS state (Eq. H14|) ~). An alter- 
nate derivation of this result is presented in Appendix 
A, where the effect of the projection operator on particle 
number fluctuations is calculated. 



Consider two BCS states defined by 
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From Eq. (fT^ . it is clear that the projection operator 
reduces the ratio of the weights of N + 2 and N particle 
states in a BCS wave function by a factor g t . Then, it 
follows that 
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The average particle number of the state I^bcs) ^ s 
given by, 



(r) 
BCS 



AM* 



(r) v 
BCS/ 



BCSl W BCS 



2 E K 



W|2 



(21) 



Since the particle numbers of |^bc S ) an d -^I^bcs) are 
identical, we can use Eq. (|18|l in Eq. I|21|l to obtain, 
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Note that g t is specified by the particle number after 
projection, A^ aftcr (= AK r )). 

Now, since the number of particles in the state I^bcs) 
before projection is given by 

^before = 2 E l^l 2 ' 
k 

Eq. 11' I'D provides us with a way to calculate the num- 
ber of particles in the state -PI^bcs) after projection, 
if I^bcs) (i-e., u k and v k ) is specified before projection. 
Eq. 122|l can be solved self consistently for ./Vafter- We 
solve Eq. (1221) numerically on a square lattice, using the 
standard BCS expressions for a d-wave superconductor, 
u k (v k ): 



where, 
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FIG. 1: (Color online) Particle density before projection 
^before aS a function of the particle density after projection 
Waiter for d-wave states with different parameters Ao- The 
dashed line indicates the Fermi liquid result ribcforc = n a fter- 



The only free parameters are the chemical potential fi 
and the variational parameter A . 

By fixing the parameter Ao, we determine the particle 
numbers (before and after projection) for various chemi- 
cal potentials. The results for particle density (n = N/L) 
are shown in Fig.^ The results clearly show that the par- 
ticle density before projection attains its maximal value 
(^before = 2), if n a f tcr = 1 (half-filling). This result holds 
for any finite value of the variational parameter Ao. In 
the opposite limit, viz., low density of electrons, ^before 
converges to the value of n a fter as expected. The size of 
the intermediate region depends on the magnitude of the 
parameter Ao, as illustrated by the results in Fig.^ 

The accuracy of Eq. IL'l'li can be checked by compar- 
ing our results with those of Yokoyama and Shiba (YS), 
who performed VMC studies of projected BCS wave func- 
tions with fluctuating particle number (but without the 
fugacity factor)i£. They determined the particle den- 
sity of the projected d-wave state P |>Fbcs) as a function 
of the chemical potential /i and the variational param- 
eter Ao, within a grand canonical scheme. The unpro- 
jected wave function I^bcs) is specified as usual, through 
Eq. (J2SJ-Eq. J57J). Since YS do not include a fugacity 
factor in their definition of the BCS wave function, pro- 
jection changes the particle number. So, we use Eq. 11' I'D 
to determine n a ftor which we compare with their results 
for particle number. 

As seen in Fig. our results are in good qualitative 
agreement with YS. Discrepancies are mostly due to finite 
corrections. YS use 6x6 and 8 x 8-lattices, while our 
analytic calculations are for the thermodynamic limi1ji&. 
The results show the singular effect of the projection near 
the insulating phase (half filling) . The chemical potential 
goes to infinity in this limit. 

In Appendix A, we present an alternate derivation of 
Eq. (1221 by a saddle point approximation without using 
a fugacity-corrected wave function. This approach also 




FIG. 2: (Color online) The particle density after projection 
n-aftcr as a function of parameter Ao for a d-wave BCS state 
at various chemical potentials fi. The figure shows a compar- 
ison between results from Eq. 12211 (solid lines) and the VMC 
results of Yokoyama and Shiba 15 (for 6x6- (circles) and 8x8 
-lattices (squares)). Numbers in the figure denote the chemi- 
cal potentials of the corresponding curves and data points. 



allows for the calculation of particle number fluctuations 
a 2 /L. 



IV. GUTZWILLER APPROXIMATION IN THE 
CANONICAL AND GRAND CANONICAL 
SCHEMES 

In this section, we discuss the differences between the 
Gutzwiller approximation in the canonical and grand 
canonical schemes. The validity of our statements can be 
checked by a comparison to nearly exact VMC— 

Let us first consider the canonical case. Here, we are 
interested in the expectation value of an operator O cal- 
culated with a particle number conserving projected wave 
function P/yP | "Fees)- The corresponding Gutzwiller ap- 
proximation can be understood as follows: 

QEbcsIPPv 6 Pn P|^bcs) 

(*BCS|PPV P|*BCS) 

(^BCSIPV 6 Pvl^BCs) 
~ 9 (*BCS|PV | *BCS) 

(^BCS|Q|^BCS) (m 
9 <*BCS|*BCS) ' 1 ' 

where Pn is the projector on the terms with particle 
number N. The Gutzwiller factor g, corresponds to the 
operator O. The first row represents a quantity which can 
be calculated exactly by fixed particle number VMGi^. 
Since the particle number is fixed, the Gutzwiller ap- 
proximation can be invoked, leading to the second row. 
The equality to the third row is guaranteed only if N is 
equal to the average particle number of | vT/bcs) (N — N). 
Here, we perform a transformation from a canonical to a 
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grand canonical ensemble, which is valid in the thermo- 
dynamic limit. 

In the grand canonical scheme, where we calculate 
the expectation value of O with a particle number non- 
conserving wave function, this scheme must be modified 
as follows: 

(^ B %\POP\^ B %) ^ (ttBCs|6|ttBCs) (2Q) 
(*BCSI PP I*BCS> " (*BCs|*BCs) ' 

where P\ ^^scs) 1S ^ ne projected d-wave state corrected 
for fugacity, i.e., a fugacity factor is included simultane- 
ously with the projection (see Sec.^J. This correction is 
essential to guarantee the validity of the Gutzwiller ap- 
proximation; without it, the lhs and rhs of Eq. I|29[) would 
correspond to states with different particle numbers. 
Comparing Eq. ffify and Eq. f25|) we get, 

(y BCS \PP N dP N P\* BCS ) ^ (^BCSI^I^BCS) 
(*BCS|^^|*BCS> ~ (*BCSI P2 |*BCS) 

(30) 

Eq. 1|29|) and Eq. (|3(J[1 constitute the main results of this 
section. Eq. (|2"9"|) shows that when the Gutzwiller approx- 
imation is used for a wave function which does not have a 
fixed particle number, a fugacity factor must be included 
along with the projection. Eq. (|30fl shows that to obtain 
identical results, one has to use different wave functions 
in the grand canonical (rhs) and canonical (lhs) schemes. 
The wave function | ^gcg) is a d-wave state corrected 
by the fugacity factor, whereas | \&bcs) is a pure d-wave 
state. Our arguments leading up to Eq. lt2$|) and Eq. (|5U|) 
can be verified by a comparison with VMC studies. We 
now proceed to do so. 

The expectation values in the canonical and grand 
canonical schemes can be calculated (nearly exactly) by 
VMC studies. In Fig. [31 we show VMC results from Gros 4 
(fixed particle number VMC, canonical) and from YS 15 
(grand canonical VMC). The discrepancy between the 
two sets of results can be explained readily by Eq. (|3U[) . 
In the case of Ao — * 0, there is only small room for 
particle number fluctuations even in the particle non- 
conserving wavefunction. Then, canonical and grand 
canonical schemes should give identical results. The 
VMC calculations in Fig.|3|do not exactly show thit be- 
havior since the grand-canonical scheme becomes inaccu- 
rate in this limit 1 ^. YS consider a pure d-wave state, i.e., 
the fugacity factor is not included in their calculations. 
In their paper, YS argued that the discrepancies between 
the two results can be removed by introducing an addi- 
tional variational parameter a, so that au =Vk/uk is 
replaced by = avk/uk (Eq. 4.1 in Ref. 15). We 
opine that the parameter a is directly related to our fu- 
gacity factor, i.e., a — gt in the wave function | ^B^g). 
This conclusion is supported by the comparison of VMC 
data to the corresponding Gutzwiller approximation (see 
below). 

The validity of the approximation in the canonical case 
(Eq. J2HJ) is weu accepted. It is used for instance, in the 



renormalized mean field theory (RMFT) of Zhang et at, 
where all physical quantities are calculated using unpro- 
jected wave functions and the corresponding Gutzwiller 
renormalization factors^, A comparison with VMC stud- 
ies with fixed particle number shows good agreement 5 
(also illustrated in Fig. 02). 

To compare the grand canonical VMC of YS with the 
Gutzwiller approximation, we need to modify Eq. H29|) . 
This is necessary because YS do not include the fugacity 
factor in their considerations, as pointed out earlier. We 
modify Eq. (|30[1 by the following procedure: 

(i) we start with a d-wave BCS state | ^bcs) for specified 
values of Ao; 

(ii) we use Eq. I|22|) to determine the chemical potential 
This fixes the particle density 

^aftcr of p|*bcs); 

(iii) we remove the fugacity factor to get | ^gcg) via Eq. 

(jTBjl . The fugacity factor is determined for n a ft er - | *bcs) 
and P\ ^bcs) correspond to the same particle density 

Rafter- 

(iv) The expectation values of the wave function 

P\ ^bcs) can now be approximated by | ^bqs) and 
Gutzwiller factors, viz., 

(*BCS|PQP|*BCS) „ (^slOl^s) 
(*BCS|PP|*BCS> ~ 5 ' 

This Gutzwiller approximation (GA) generalizes Eq.[21 
for wave functions that do not conserve particle number. 
In Appendix B, we discuss this approximation for the 
different terms in the t — J model. 

In Fig. 0] we compare the GA of the kinetic energy 
and the expectation value E^ 2 \ of the remaining 
terms in the t — J model ((SiSj), (hifij), and the 3-site 
term) to those from the grand canonical VMC. A good 
agreement between the VMC and Gutzwiller results is 
seen, which confirms the validity of our grand canonical 
Gutzwiller approximation (Eq. H29f lV 

In Fig. Inland Fig. 21 we also show Gutzwiller approx- 
imations for the fixed particle number VMGi Clearly, 
canonical and grand canonical approaches yield different 
energies (as do the corresponding VMC studies ). Wc 
emphasize this is because of the projection operator P, 
which changes the particle number in a grand canonical 
scheme. For these two methods to yield the same re- 
sults, a fugacity corrected wave function must be used 
when working in a grand canonical ensemble. Hence, all 
previous speculations about the coincidence of these two 
VMC schemes in the thermodynamic limit have to be 
reformulated carefully. 

V. SUMMARY 

In this paper, we considered the effects of Gutzwiller 
projection on a state which does not have fixed parti- 
cle number. We showed that it is necessary to include a 
fugacity factor when invoking the Gutzwiller approxima- 
tion for such states. The effects of projecting a number 
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FIG. 3: (Color online) (a) The kinetic energy E {1) and (b) 
the energy of the remaining terms E^ 2 ' per site of the t — 
J model as a function of the variational parameter Ao for 
the d-wave state at a filling n = 0.9. Fixed particle (can.) 
VMC data- (circles, 82 sites) and grand canonical VMO— 
(squares, 8x8 sites) are compared. The dashed/solid lines 
represent the corresponding Gutzwiller approximations (GA). 
For a detailed description, we refer to the text and Appendix 
B. 




FIG. 4: (Color online) (a) The kinetic energy E w and 
(b) the energy of the remaining terms E^ per site of the 
t — J model, as a function of the variational parameter Ao 
for the d-wave state at various densities. The data points 
(diamond, circle, square) are taken from the grand canonical 
VMC study of Yokoyama and Shiba»». The solid/dashed lines 
represent the corresponding Gutzwiller approximations for a 
grand canonical/fixed-particle VMC study. For a detailed de- 
scription, we refer to the text and Appendix B. 



non-conserving BCS state were studied by examining the 
relation between particle number before and after projec- 
tion. We obtained an analytical expression, Eq. (|22p. and 
compared to variational Monte Carlo data (Fig. [SJ. We 
discussed the discrepancies in the VMC results for pro- 
jected BCS wave functions obtained in the canonical and 
grand canonical schemes, and presented a resolution. In 
conclusion, we have clarified several subtle properties of 
the Gutzwiller projection operator P acting on a BCS 
state, and hope that these results lead to a better under- 
standing of the Gutzwiller approximation in the grand 
canonical scheme. 
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by the Deutsche Forschungsgemeinschaft. V. N. M. ac- 
knowledges partial financial support from The City Uni- 
versity of New York, PSC-CUNY Research Award Pro- 
gram. 



APPENDIX A: SADDLE POINT 
APPROXIMATION TO MEAN PARTICLE 
NUMBER AND NUMBER FLUCTUATIONS IN 
PROJECTED BCS WAVE FUNCTIONS 

In Sec. Ill, we used the fugacity factor to derive 
Eq. 1)22(1. Here, we present an alternative approach by a 
saddle point approximation to discuss the effects of pro- 
jection on the mean particle number of a BCS state. This 
approach also describes the particle number fluctuations 
after projection. 

The particle number distribution for an unprojcctcd 



BCS wave function 



(o) 
Pn = 



(0) 



Pn 



can be written as 



dA 



m(£) no 



M 2 A) 
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where N a — N/2 is the number of electron pairs. This 
relation can be checked by expanding the product in the 
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wave function 

|4-BCS 



> =n( Mfe+wfec iT c -fc,i) i°) ' 

k 



and considering contribution to (^bcs I^bcs) from each 
term. In Sec. Ill, we showed that the particle number 
distribution of a projected wave function is related to 
the unprojected distribution pff by 



(o) 

Pn = 9n P N , 



(A2) 



where, 



9n ~ C 



{{L-N/2)\f 
L\{L-N)\ 



Particle number and number fluctuations of the projected 
wave function can be derived from Eq. I|A2(1 and Eq. ilAlfl . 
upon invoking a saddle point approximation. We define 
a generating function 5a, 



2n(KI 2 + kl 2A ) 



by 
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d\ dX X 
_ d_ 

~ dX 



Ft — -I- 1 

^log(| Ufc | 2 + k| 2 A)-^ + - 



E 



kl 2 f + i 



u k \2 + \v k \ 2 X X 



(A9) 



for N > 1, 



n w 2 ^ 

L 



AM 



L V kl 2 + AM 2 



(A10) 



and 



<9f 1 TV 

^ = - - log A - log(L - -) + log(L - N) 



, 



(All) 



E 



(A3) 



We invert Eq. <| A3|> using a contour integral on the com- 
plex A-plane along a circle around A = 0: 



(0) 



Pn 



>— 'A 



1 



2tTC / A JV » + 



-dA . 



(A4) 



Note that, in the integrand, only p^/v M = Pn gi yes 
a finite value. The others powers of A vanishes. Multi- 
plying by g N gives 



PN^^-JgN^dX 
Eq. i|A5() can be written as 

p N ^J-ld\ e /(W) 

27TI 



(A5) 



(A6) 



where 



,N 



f{X,N) =]ogS x -(— + l)logA +log ffw . (A7) 

Using Stirling's formula, 

TV iV 
loggjv ~ 2 (L — — ) l g(£--)-(i-iV) log(L-iV) 



-L logL +logC 



(A8) 



The saddle point (n = A) of Eq. I|A6|I is determined 



L-N 



N 



=9t- 



(A12) 



Eq. ifATOjl and Eq. llAT^ lead to Eq. (j^ . 



,2|„„ 12 



- ^ 1 ^ g \v k 

Rafter = n = — = 2 



for the average particle density of a projected BCS wave 
function. Without the factor g^r this calculation would 
give the well known result for an unprojected BCS wave 
function, 

"before = n= 2 — ^2 k| 2 ■ 

k 

To calculate the particle number fluctuations of the 
projected wave function, we need to expand f(X,N) up 
to second order in N and A around the saddle point. 
Then, integration over A in Eq. I)A6|I approximates the 
particle number distribution p by a gaussian distribution, 
yielding an expression for number fluctuations. With 



Vk 



N 



fx N (X,N) 
/nn(X, N) 



V (KI 2 + kl 2 A) 2 A 2 

d 2 f 1_ 

dXdN 2A ' 

d 2 f _ 1 1 1 

diV 5 ~ 2L- ? ~ L — N ' 
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FIG. 5: (Color online) The fluctuation a^ ftcI / o"b c f orc af- 
ter/before projection for different values of the variational 
parameter Ao = 0.01/0.1/1/10, as a function of the particle 
density after projection n a ft C r. 



the second order expansion can be written as 

(A -A) 2 



/(A, A) -/(A, A) «/aa(A,A0- 



+/aw(A,A)(A-A)(A- 
+f NN CX,N) {N ~ N)2 



A) 



(A13) 



For this level of the saddle point approximation for A, 
the contour around A for the integral in Eq. (|A6ll must be 
taken so that / AA (A, A) (A- A) 2 < 0. Since f xx (A, A) > 
and contribution only near the saddle point is relevant, 
the path is taken from A — ioo to A + ioo. By variable 
transformation A = A + iX' , one can perform a gaussian 
integral of A'. Then, we obtain a a gaussian distribution 
for pn, 



Pn 



1 

2~~ 



— / dX' e /(*+^) 



exp 



NN 




3 /(A,iV) 



(A14) 



The variance of A (average particle number of a pro- 
jected BCS wave function) can now be read off from 



Eq. ||A14|) . We get, 



J afte 

L 



Inn 



= 2 



1 



(l-g)(l-n) 



^l"fcl 2 l^fc 



l l^k (KP+Kp A) 2 



(A15) 



where we used Eq. (|A10|) in the second term. For A we 
must insert g\. For completeness, we mention that for 
the unprojected wave function, i.e. not included, this 
approach yields the known result 



before 



= 4 



Ufc| 2 |Vfe| 2 



(A16) 



2 __ 

The fluctuations ^- are illustrated in Fig.|5]as a function 
of the particle density after projection n a fter (= n) for 
unprojected (| ^bcs)) an d projected (P\ *bcs)) BCS d- 
wave functions. As expected, the fluctuations vanish at 
half filling, since projection freezes the charge degrees of 
freedom entirely. 



APPENDIX B: GUTZWILLER 
APPROXIMATION FOR THE t J 
HAMILTONIAN 

We summarize the Gutzwiller approximation for the 
so-called 3-site terms in the t—J model, that are included 
in the VMC study of Yokoyama and Shiba— . 

The t—J model can be derived from a large U expan- 
sion of the Hubbard model. The Hamiltonian is valid in 
the reduced Hilbert space of no double occupied states, 
and is given by 



fleff = T + H 



(2) 



eft 



where 



T = -t E ( c L c J> + c L Ci 



(Bl) 
(B2) 



and 



H, 



(2) 



4E 



{i,3) 



J 



E 

E 4 



(B3) 



t 2 

Here, J = 4jj, Si are the spin operators on site i, and 



( f + rii , with rii tC 



are pairs of n.n 



sites and i + r denotes a n.n. site of i. 
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We are interested in the energies and E<-V calcu- 
lated in Ref. Ill 



L (Vbcs\P P\^bcs} ' 

1 (fgciPjl P^BCs) 
i (*BCS|-P .PI^BCS) 



(B4) 



We invoke the Gutzwillcr approximation. The renormal- 
ization factors, gt for kinetic energy (Eq. i|B2l) ) and <?s 
for spin exchange (first term in Eq. (|B3(l ). are given in 
Eq. ©. The second term of Eq. (|B3|) . riiUj, is not renor- 
malized. The approximation for the 3-site terms (3rd and 
4th term of Eq. (JB3J)) is done as follows (\<ip) = P\ipo)): 

(V\4 +T ,icl l a, l c i+Tl ,- [ \V) 
(*l*> 

(*|*) 

(*o|c|+r T ni,4.(l - n l!T )c i+r ,, T |* ) 
= 53 LJ /lTr ,, T , x , (B5) 



(*o|*o) 



(y\cl +Td 4 A c ld c l+T ,,^) 



93 



(*o|*o) 



(B6) 



The renormalization factor g% is derived by considering 
the number of terms that contribute to the projected 
and the unprojected side respectively. The projected 
side (lhs) contributes only if (i) site i + r is unoccu- 
pied, i.e., probability (1 — n), (ii) site i is singly occu- 
pied by a ^.-electron, i.e., probability ni, and (iii) site 
i + t' is singly occupied by an f -electron, i.e., probabil- 
ity 7ij. On the other hand, the unprojected side (rhs) 
in Eq. JBSJ/Eq. I|B6() contributes only if (i) site i + r 
is not occupied by an |-electron/|-electron, i.e., proba- 
bility (1 — Tif) / (1 — nj.) , (ii) site i is singly occupied 
by a J,-electron, i.e., probability n|(l — n-j-), and (iii) site 
z + t' must have an "("-electron, i.e., probability n-\. These 
probabilities yield the Gutzwiller factor (ratio of contri- 
butions from projected and unprojected states), 



93 



(1 - n)n a n a 



1-n 

(1 - n a )n a {l - n a )n a (1 - n a ) 2 



(B7) 



where we assumed raj = n± = n CT . 

We can now write down the renormalized £ — J Hamil- 
tonian H' P 



off ' 



(B8) 



where, 



t' = -g t t Yl (4^>+ c L 



(B9) 



(<,i> 



-53 



+93 J C \+T,-A,* C t,-V C t+T'.<7 • (BIO) 



By using V and Eq. (JB4J) (£jW and £( 2 )) can be 

calculated using unprojected wave functions, 



e{1) = 1 (^BCSIT'I^BCS) 

L (*BCs|*BCS> ' 

E(2) _ 1 (*BCS|^|*BCS) 



L *BCS * 



BCSl^BCS/ 



(Bll) 



Evaluating Eq. (|B11|) by Wick's decomposition for a 
d-wavc BCS state, 



-2.9t(^+^) 



ff( 2 ) 
0.25J 



(B12) 



(3 5s -l)(£ 2 +£ 2 )/2 

-(3fo + l)(|A s | 2 + |A B | a )/2 - 2n 2 

- 93 n (1 - n ff ) (6a + 6w + 2£e-i/ + 2 ^+ a ) 
-53(1 + ^)^+^+4^) 

- ff3 (2 - n„) (IA.I 2 + |A,| 2 + 45R(A ;c A;)) , 
where we defined (r, r' = x, y) 

Cr = ^(cj^Ci+^o.) = i^2cOs(fc r )|w fc | 2 , 
a k 

£r±r' = Y(4+t,<j C i±t' ,a) = ^X2cOs(fc T ±/c T 0|Wfc| 2 , 
cr k 

A T = (cJ ;T cJ +t4 - ct 4 ct +TiT ) = ~^2cos(fc T H"ft . 



The last three rows in Eq. (|B12|) correspond to the 3-site 
terms of the t — J model and are renormalized by the 
Gutzwiller factor 173. Here, it's important to note that 
the order parameter A r is related, but not identical, to 
the previous introduced variational parameter A . 
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